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Abstract 


The  paper  outlines  existing  knowledge  on  A-optimal  designs  for  comparing  test  treat¬ 
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1.  Introduction 

— We-wottki^likeTo  introduce  the  problem  with  an  example.  How  should  we  design  an 
experiment  to  compare  4  test  treatments  with  a  control,  using  18  experimental  units?  As 
a  statistical  question  we  will  not  be  able  to  answer  it  unless  it  is  asked  in  a  more  precise 
manner.  To  begin  with  we  need  to  postulate  a  model  for  the  response  observed  upon 
application  of  a  treatment,  test  treatment  or  control,  to  an  experimental  unit.  '  In  this 
papery we-shaU  consider  three  possible  models: 

O-way  elimination  of  heterogeneity  model  in  which  all  experimental  units  are  homogeneous 
before  application  of  treatments^ 


Vxj  =  M  +  n  +  ei} 


(1.1) 


1-way  elimination  of  heterogeneity  model  in  which  the  experimental  units  can  be  divided 
into  several  homogeneous  blocks: 


Vij  =  H  +  U+  Pj  +  tij 


(1.2) 


t  2-way  elimination  of  heterogeneity  model  in  which  the  experimental  units  can  be  concep¬ 
tually  arranged  according  to  rows  and  columns: 


Vijl  =  H  +  U  +  pj  +  pt  +  Ujt  (1.3) 

Now  we  can  be  more  precise  about  what  we  mean  by  comparing  test  treatments  with 
a  control.  Our  goal  is  to  estimate  the  magnitude  of  each  (tv  —  To).  Assuming  that  the 
error  component  e  in  the  model  is  homoscedastic,  the  method  of  least  squared  will  be  used 
to  estimate  the  contrast  (r,  —  r0)  ;  this  happens  to  be  the  best  linear  unbiased  estimator 
(f,  —  ?o)  •  In  assigning  the  treatments  to  the  experimental  units  we  have  to  make  sure  that 
the  contrasts  (t,  -  To)  are  estimable.  In  case  we  have  more  than  one  choice  for  making 
this  assignment  we  want  to  select  the  one  which  guarantees  high  efficiency  in  the  sense  of 
achitving  the  minimum  value  of 
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YlVar 

«  =  1 

A  design  which  gives  the  minimum  will  be  called  an  A-optimal  design. 

Without  further  ado  we  give  A-optimal  designs  under  each  of  the  three  models. 
A-optimal  design  under  model  (1.1): 

Assign  3  experimental  units  to  each  of  4  test  treatments  and  6  to  the  control. 

A-optimal  design  under  model  (1.2),  when  there  are  6  blocks  of  size  3  each: 

Take  each  column  of  the  following  array  as  a  block: 


0  0  0  0  0  0 
1112  2  3 
2  3  4  3  4  4 


Here  0  denotes  the  control  and  1,2, 3,4  the  test  treatments  . 


A-optimal  design  under  model  (1*3),  when  there  are  3  rows  and  6  columns: 
Assign  the  treatments  according  to  the  following  array: 


1  0  3  4  2  0 
0  3  4  2  0  1 
4  2  0  0  1  3 


Even  this  small  size  example  demonstrates  the  fact  that  the  problem  of  finding  an  A- 
optimal  design  is  a  difficult  one.  During  the  past  several  years  there  has  been  a  concerted 
effort  to  identify  and  construct  A-optimal  designs  for  the  general  problem  of  comparing  v 
test  treatments  with  a  control.  In  this  paper  we  shall  attempt  to  summarize  these  results 
which  we  hope  will  be  useful  to  both  the  theoretician  and  the  practitioner. 

In  section  2  we  give  general  results  for  A-optimal  designs  for  comparing  v  test  treat¬ 
ments  with  a  control  in  each  of  the  3  models  (l.l),  (1.2)  and  (1.3).  In  section  3  we  give 
model  robust  A-optimal  designs.  In  section  4  we  suggest  various  approaches  for  finding 
efficient  designs  in  those  cases  where  A-optimal  designs  are  unknown.  In  section  5  we  give 
A-optimal  designs  for  comparing  test  treatments  with  two  or  more  controls.  In  section 


6  we  give  an  overview  of  the  literature  of  optimal  designs  for  comparing  test  treatments 
with  controls. 


2.  A-optimal  Designs. 

We  shall  give  A-optimal  designs  for  comparing  t;  test  treatments  with  a  control 
separately  for  the  O-way,  1-way,  and  2-way  elimination  of  heterogeneity  model.  Throughout 
this  section  the  control  will  be  denoted  by  the  symbol  0  and  the  test  treatments  by 
1, 2,  ...,v  . 

2.0.  A-optimal  designs  for  0-way  elimination  of  heterogeneity. 


Our  statistical  setup  consists  of  n  experimental  units,  and  our  model  of  response  under 
a  design  d  is: 


y d ij  —  P  "b  Ti  +  (2-0) 

where  j  •=  1,...,^,,  t  =  0,1,..., v  .  Here  rdi  is  the  number  of  experimental  units  receiving 
tieatment  i  .  We  assume  the  model  to  be  additive  and  homoscedastic  .  The  symbols  in 
equation  (2.0)  have  their  standard  meaning.  The  A-optimal  design  minimizes 

V 

5^(l/rd0  +  l/rdi) 
i=  1 

subject  to  the  restriction  rd0  +  rdi  + . . .  +  rdv  —  n  .  In  case  v  is  a  square  and  n  =  m  (u  +  v^) 
for  an  integer  m,  the  A-optimal  design  d*  is: 

r<j*i  =  •  *  •  =  rj‘V  =  m,  rd~o  =  my/v. 

2.1.  A-optimal  designs  for  1-way  elimination  of  heterogeneity. 

Our  statistical  setup  consists  of  6  blocks  of  size  k  each,  and  the  model  of  response 
under  a  design  d  is: 


Vdijp  —  M  “b  "b  0}  ~h  fyp  (2T) 

where  *  =  0, 1, . . . ,  v  ,  j  =  1, . . . ,  b  and  p  =  0, 1, . . . ,  .  Here  is  the  number  of  times 

treatment  i  is  used  in  block  j.  Let  Nd  denote  the  matrix  (n^.y  )  , 


Cd  =  Diag(rd0,rdu...,rdv)  -  k  lNdN'd, 


P  being  atxv  +  l  matrix.  Then  an  A-optimal  design  minimizes 

trace  P  CJ  P'  (2.2) 

over  all  possible  designs. 

Experience  has  shown  that  this  minimization  is  a  difficult  task.  As  in  other  cases  of 
exact  design  theory,  it  is  highly  unlikely  that  we  can  obtain  one  method  which  is  capable 
of  producing  A-optimal  designs  for  arbitrary  values  of  v,b  and  k  .  Recently  several  families 
of  A-optimal  designs  have  been  discovered. 

At  this  point  it  is  useful  to  recall  a  celebrated  result.  If  there  was  no  control  and  if  we 
were  interested  in  comparing  v  test  treatments  among  themselves  then  a  BIB  design  in 
the  v  test  treatments  would  be  A-optimal.  Unfortunately,  with  the  presence  of  the  control 
and  for  the  set  of  contrasts  of  interest  a  BIB  design  is  almost  never  an  A-optimal  design. 
However,  we  can  sometimes  utilize  BIB  designs  in  the  test  treatments  to  construct  an 
A-optimal  design  for  our  problem.  We  shall  give  some  families  of  such  designs  below.  For 
convenience,  we  introduce  the  notation  ABIB  (v,6,  k  -  t\t)  to  denote  a  BIB  design  in  the 
v  test  treatments  in  b  blocks  of  size  k  -  t  each  augmented  by  t  replications  of  the  control 
in  each  block. 

Family  1.  An  ABIB(v,6,fc  -  1,  l)  is  A-optimal  whenever  (k  -  2)2  +  1  <  v  <  (k  —  l)2  . 

An  example  of  an  A-optimal  design  when  v  =  7,6=7,k  =  4is  given  below,  where 
the  columns  are  the  blocks: 


0  0  0  0  0  0  0 

1  2  3  4  5  6  7 

2  3  4  5  6  7  1 

4  5  6  7  1  2  3 

For  each  (v,Jfc)  satisfying  ( k  —  2)2  +  1  <  v  <  (k  -  l)2  ,  there  are  an  infinite  number 
of  A-optimal  ABIB(v,i>,  k  —  1;  l)  designs.  These  results  and  more  details  are  available  in 
Hedayat  and  Majumdar  (1985a). 

Stufken  (1986)  has  generalized  the  preceeding  idea  to: 

Family  2.  An  ABIB(v,  b,k-t\ t)  is  A-optimal  whenever  (A:  —  t  -  l)2  +  1  <  t2v  <  (k  -  t)2  . 
An  example  of  an  A-optimal  design  when  v  =  8,  b  =  28  ,  k  —  8  is  given  below: 
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Sometimes  we  can  use  two  BIB  designs  to  construct  an  A-optimal  design  for  our 
problem.  We  give  below  one  such  family,  which  is  taken  from  Cheng,  Majumdar,  Stufken 
and  Ture  (1986): 

Family  3.  For  v  =  a2  —  1,  6  =  7(a  +  2)(a2  -  1)  and  k  —  a  ,  the  union  of  an  ABIB(t;,7(a-)- 
l)(a2  -l),a  -  1;  1)  and  a  BIB  design  in  all  the  v  +  1  treatments,  test  treatments  and 
control,  in  ~ia(a  -f  1)  blocks  of  size  k  each  is  A-optimal  whenever  a  is  a  prime  power,  and 
7  is  any  integer. 

An  example  when  a  =  3,7  =  l,v  =  8,6  =  40,fc  =  3is: 


;• "OP 


00000000000000000000 

11111112222223333344 

23456783456784567856 


00000000147123123123 

44555667258456564645 

78678788360780078807 

Stufken  (1986)  has  some  more  families  of  A-optimal  designs. 

To  establish  the  optimality  of  these  families,  the  starting  point  is  a  result  due  to 
Majumdar  and  Notz  (1983),  which  we  now  proceed  to  state.  There  are  many  parameter 
combinations  (v,b,k)  which  do  not  belong  to  any  of  the  three  families  but  for  which 
Majumdar  and  Notz’s  result  could  still  be  used  to  get  an  optimal  design.  A  complete  list 
of  all  designs  available  by  this  result,  when  2<A:<8,A:<v<30,v<6<50is  given  in 
Hedayat  and  Majumdar  (1984).  Before  stating  the  result  we  need  some  definitions. 

DEFINITION  2.1.  d  is  a  Balanced  Treatment  Incomplete  Block  (BTIB)  design  if 

A d01  =  • • •  =  ^dOv 
A(ji2  =  •  =  ^dv-lv 

where  A dij  —  5Zp=i  ndipndjp  ■  This  definition  is  due  to  Bechhofer  and  Tamhane  (1981). 

DEFINITION  2.2.  For  integers  t  G  {0,1,...  k  -  1}  and  s  G  (0, 1, . . .  ,b  —  1}  d  is  a 
BTIB[v,b,k]t,s)  if  it  is  a  BTIB  design  with  the  additional  property  that 

Tldij  G  {0,  l} ,  1  1,  .  .  .  ,  V,  j  1, . . . ,  6. 

rtrfOl  =  •  *  *  =  KdOs  —  t  +  1 
**d0a+l  =  ■  •  ■  =  Tldob  =  t. 

A  BTIB(v,b,k;t,s)  is  called  a  Rectangular  (R-)  type  design  when  s  =  0,  and  a  Step  (S-) 
type  design  when  s  >  0.  The  layout  of  these  designs  can  be  pictured  as  follows,  with 
columns  as  blocks,  in  each  of  the  two  cases  R-type  and  S-type: 
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di  and  d 2  are  components  of  the  design  which  involve  the  test  treatments  only. 


He  shall  also  need  the  following  notations: 


a  =  (v  —  l)2,  c  ~  bvk(k  —  1),  p  =  v(k  —  1)  +  k, 

A  =  {(x,  z)  :  x  =  0, . . . ,  [k/2]  —  1;  z  =  0, 1, . . .  ,6  with  z>  0  when  x  =  0}. 
Here  [•]  denotes  the  largest  integer  function. 


g( x,  z)  —  a/{c  -  p{bx  +  2)  +  (6x2  +  2x2  +  2)} 

+  \/{k(bx  +  2)  -  (6x2  +  2x2  +  z)}. 

Now  we  are  ready  to  state  the  result  of  Majumdar  and  Notz(l983): 

THEOREM  2. 1 .  Let  v,b,k  be  integers  with  k  <  v  .  A  BTIB(v,b,k\t,s )  is  A-optimal  in 
the  class  of  all  designs  if 


g[t,s)  =  Min{g{x,z)  :  [x,z)  G  A}. 

Hedayat  and  Majumdar  (1984)  have  devised  an  algorithm  for  obtaining  A~optimal 
designs  based  on  Theorem  2.1.  Jacroux  (1986)  has  generalized  this  algorithm.  His  algorithm 
is  often  capable  of  producing  A-optimal  designs  which  are  not  necessarily  BT1B  in  their 
structure.  An  example  is  a  group  divisible  partially  balanced  block  design  in  the  test 
treatments  ,  extended  by  one  replication  of  the  control  in  each  block.  Below,  we  give  one 
such  example  of  an  A-optimal  design  when  v  =  9,  6  =  9  and  k  =  4: 

000000000 
1  1  1  2  2  2  3  3  3 

456456456 
789897978 

2.2.  A-optimal  designs  for  2-way  elimination  of  heterogeneity.  Our  statistical  setup  con¬ 
sists  of  bk  experimental  units  arranged  in  a  k  x  b  array,  and  the  model  of  response  under 
a  design  d  is: 


(»'  =  0,1,..., v;  j 
unit  in  cell  (£,/). 


Vdijt  ~  P  +  U  +  Pj  +  Pl  +  t\jt  (2.3) 

=  l  =  1  if  treatment  i  is  applied  to  the  experimental 


ndij  =  number  of  times  treatment  *  occurs  in  column  j, 

tUdit  =  number  of  times  treatment  i  occurs  in  row  i , 
b 

rdi  =  ^2  ndij 

Nd  =  (Mrfiy) ,  a  v  x  b  matrix 

Md  =  (m<it/) ,  a  v  x  A:  matrix 

P  is  the  u  x  (v  +  l)matrix  defined  in  subsection  2.1. 

f’d  =  (rdChrdl  j  •  •  •  i  •’fiu) 

Cd(2)  =  (rdo,rdi,...,rdv)  -k^NaN^-  b~lMdM'd  +  (fc*)"1^ 

Then  an  A-optimal  design  minimizes  trace  PC^^P'  .  We  shall  now  highlight  some  of  the 
results  from  recent  literature. 

Family  1.  Let  p  be  an  integer  and  v  =  p2,b  =  k  =  p2+p.  A  b  x  b  array  in  which  each 
test  treatment  appears  once  in  each  row  and  in  each  column  and  the  control  appears  p 
times  in  each  row  and  in  each  column  is  A-optimal. 

One  easy  way  to  construct  members  of  this  family  is  to  start  with  a  Latin  square  of 
order  p2  +  p  and  change  symbols  p2  +  1 , . . . ,  p2  +  p  to  0  (control) .  We  illustrate  this  in  the 
following  example  with  v  =  4,b  —  k~6. 
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This,  and  some  more  general  results  are  available  in  Notz  (1985). 
This  has  been  generalized  by  Majumdar  (1986). 


Family  2.  Let  p,  a  and  7  be  integers,  v  —  p2  ,  k  =  a(p2  +  pi  and  b  =  nf(p2  +  P)  •  A 
k  x  b  array  in  which  each  test  treatment  appears  a  times  in  each  column  and  7  times  in 
each  row,  and  the  control  appears  ap  times  in  each  column  and  7p  times  in  each  row  is 


A-optimal. 


One  way  to  construct  members  of  this  family  is  to  form  the  array 


(L,y),  *  =  l,...,a;j  =  1....7 


where  each  L,y  is  a  member  of  Family  1. 


Family  3.  A  k  x  b  array  is  A-optimal  if 

(i)  it  is  an  A-optimal  block  design  for  one  way  elimination  of  heterogeneity  with  columns 
as  blocks,  and 

(ii)  the  total  number  of  replication  for  each  treatment,  test  treatment  or  control,  is  divided 
equally  among  the  k  rows. 

This  has  been  given  by  Jacroux  (1984c).  The  following  is  an  example  when  v  =  9, 6  = 
24,  A:  =  3  : 


041025783096008079162345 

105802005204466907213879 

310140220730950680976458 


3.  Model  Robust  A-optimal  Designs. 

There  are  circumstances  in  which  the  experimenter  is  not  sure  whether  to  fit  a  1-way 
or  a  2-way  elimination  of  heterogeneity  model  to  the  data.  In  such  situations  it  would 
be  highly  desirable  to  obtain  a  design  which  is  A-optimal  under  each  of  these  models. 
Hedayat  and  Majumdar  (1985b)  studied  this  aspect  of  the  problem  and  gave  some  families 
of  model  robust  A-optimal  designs.  The  families  were  constructed  using  the  Euclidean 
plane,  the  Projective  plane  and  some  other  geometrical  structures.  The  exact  description 
of  the  families  are  somewhat  involved;  some  typical  examples  are  given  below. 

Example  1.  let  v  =  4  ,  k  =  3  and  5  =  6.  The  following  design  is  A-optimal  for  both 
one-  and  two-way  elimination  of  heterogeneity  models: 

1  0  3  4  2  0 

0  3  4  2  0  1 

4  2  0  0  1  3 

In  fact,  this  design  is  A-optimal  for  0-way  elimination  of  heterogeneity  model,  too. 

Example  2.  Let  v  =  7  ,  k  —  4  and  5  =  28.  The  following  design  is  A-optimal  for  both 
one  and  two-way  elimination  of  heterogeneity  models: 


0000000  1234567  1234567  1234567 
1234567  0000000  2345671  2345671 
2345671  2345671  0000000  4567123 
4567123  4567123  4567123  0000000 


Before  closing  this  section  we  would  like  to  mention  that  the  designs  in  families  1  and 
2  in  subsection  2.2  are  A-optimal  under  0-way,  1-way  and  2-way  elimination  of  heterogene¬ 
ity  models,  while  the  designs  in  family  3  are  A-optimal  at  least  under  1-way  and  2-way 
elimination  of  heterogeneity  models. 

4.  Other  Efficient  Designs. 

Even  though,  for  each  set  of  v  test  treatments  there  is  an  A-optimal  design  for  a  0-, 
1-  or  2-way  elimination  of  heterogeneity  model,  the  task  of  finding  this  design  can  be  very 
difficult  indeed.  For  situations  where  an  A-optimal  design  is  unknown,  there  are  several 
alternative  ways  of  planning  an  experiment.  Here  are  some  possiblities: 


4.1.  Modify  the  optimality  criterion,  so  that  for  the  given  circumstances,  an  optimal  de¬ 
sign  with  respect  to  the  new  criterion  can  be  constructed.  One  reasonable  alternative 
to  A-optimality  is  MV-optimality.  An  MV-optimal  design  is  the  one  which  achjetves  the 
minimum  value  of 

Maxi<i<vVar  [fdi  -  fd 0) 

among  all  competing  designs. 

For  example  when  v  =  6,  b  =  11  and  k  =  3  ,  as  well  as  when  u  =  6  ,  6  =  16  and 
k  =  4  ,  an  A-optimal  block  design  is  unknown.  However,  we  are  able  to  give  MV-optimal 
block  designs  for  these  parameters.  These  are  exhibited  below: 

Example  1.  v  =  6,  b  =  11  and  k  =  3. 

00000000014 

11122233325 

45645645636 

Example  2.  v  =  6,  b  =  16  and  k  =  4. 

0000000000000000 

1123131212121312 

2544254335432434 

3665466646655665 

It  is  interesting  to  note  that  the  design  in  example  1  is  an  S-type  BTIB  design,  while 

the  test  treatments  in  the  design  in  example  2  form  a  group  divisible  design.  These  and 

more  MV-optimal  designs  can  be  found  in  Jacroux  (1984b, 1984c  and  1985). 

4.2.  Limit  the  class  of  competing  designs  to  treasonably  rich”  subclass,  so  that  an  A- 
optimal  design  within  this  subclass  can  be  constructed.  For  example,  under  a  1-way 
elimination  heterogeneity  model  the  BTIB  designs  fctm  such  a  subclass.  For  v  =  3  , 
b  =  15  ,  k  =  2  the  A-optimal  design  is  not  available  in  the  literature,  but  a  design  which 
is  A-optimal  among  all  BTIB  designs  is  given  by 


This  approach  has  been  studied  in  Hedayat  and  Majumdar  (1984),  where  some  series 
of  such  designs  have  been  cataloged.  Numerical  evidence  indicates  that  optimal  designs 
obtained  in  this  fashion  are  highly  efficient  in  the  class  of  all  designs.  There  are,  however, 
isolated  instances  where  they  perform  poorly.  A  similar  study  for  the  MV-optimality 
criterion  has  been  carried  out  by  Jacroux  (1985). 

4.3.  Search  for  an  approximately  A-optimal  design.  This  can  be  carried  out  in  two  steps. 
First  compute 

kvg(t,s)  =  kvMin{g(x,z)  :  (i,g)  6  A}  (4.1) 

where  the  function  g(x,  z)  has  been  defined  in  section  2.1  just  before  Theorem  2.1.  This 
gives  a  lower  bound  to  the  value  of  the  A-criterion  given  in  expression  (2.2).  The  second 
step  consists  of  guessing  a  good  design  based  on  available  theory.  Compute  the  corre¬ 
sponding  value  of  the  expression  (2.2)  for  this  design  and  compare  with  the  minimum 
value  given  in  (4.1).  If  the  comparison  is  poor  in  the  opinion  of  the  experimenter,  then  he 
should  modify  his  guess  and  try  again. 

Let  us  demonstrate  this  approach  by  an  example.  Let  v  —  21  ,  6  =  30  and  k  =  9  . 
Here  the  minimum  given  by  (4.1)  is  2.589.  Our  experience  shows  that  BIB  designs  in 
the  test  treatments  augmented  by  one  or  more  replications  of  the  control  in  each  block  is 
often  highly  efficient,  as  seen  in  Families  1  and  2  of  subsection  2.1.  In  our  case  we  can 
try  a  BIB(21,30, 10,7,3)  in  the  test  treatments  augmented  by  2  replications  of  the  control 
in  each  block.  For  this  design  the  value  of  (2.2)  is  2.618,  giving  an  efficiency  of  at  least 
98.87%.  So  this  is  indeed  a  highly  efficient  design.  This  approach  of  approximating  an 
A-optimal  design  by  an  augmented  BIB  design  has  been  studied  by  Stufken  (1986). 

Another  method  of  tracking  down  a  good  approximation  has  been  given  in  Cheng, 
Majumdar,  Stufken  and  Ture  (1986).  It  consists  of  first  determining  the  point  (t,s)  which 
minimizes  the  function  g(x ,  z)  given  in  (4.1).  In  case  a  BTIB(v,  b ,  k ;  t,  s)  exists,  it  is  optimal 
by  Theorem  2.1.  If  it  does  not,  then  at  least  one  of  the  following  two  designs  is  expected 
to  be  a  good  approximation: 

(i)  A  design  with  the  same  number  (bt+s)  of  replications  of  the  control  as  a  BTIB(v,  b ,  it;  t,s) 


l.'  v 


and  “combinatorially  close”  to  a  BTIB  design. 

(ii)  A  BTIB  design  with  the  number  of  replications  of  the  control  “close”  to  bt  +  s  . 

We  demonstrate  the  idea  by  an  example  when  v  =  5  ,  b  =  7  and  k  =  4  .  Here 
(t,s)  =  (1,0)  and  vkg(t,s)  =  2.04  .  There  is  no  BTIB(5, 7,4;  1,0)  .  Consider  the  following 
two  designs: 

0000000  0000001 

,1111122  ,  0011222 
1  '  2  2  3  3  4  3  3  ’  2:1334343 

4545545  2455554 

Here  bt  +  s  =  7  ,  d\  is  a  non-BTIB  design  with  7  replications  of  the  control,  while  d2  is 
a  BTIB  design  with  8  replications  of  the  control.  The  value  for  expression  (2.2)  for  d\  is 
2.058  giving  it  an  efficiency  of  99.2%.  The  value  for  expression  (2.2)  for  d2  is  2.143  giving  it 
an  efficiency  of  95.2%.  Thus  both  of  these  designs  are  highly  efficient,  d\  being  the  better 
of  the  two. 

Finally,  with  the  availability  of  today’s  high  speed  computers  and  supercomputers, 
one  can  find  the  A-optimal  design  by  a  complete  search  among  all  possible  designs  if  the 
parameters  v,b  and  k  are  not  too  large. 

5.  A-optimal  Designs  For  Two  or  More  Controls. 

So  far  we  have  been  discussing  optimal  designs  for  comparing  test  treatments  with 
one  control.  There  axe  circumstances  when  the  test  treatments  have  to  be  compared  with 
two  or  more  controls.  Suppose  we  denote  by  S  the  set  of  all  controls  and  by  T  the  set  of 
all  test  treatments.  Then,  an  A-optimal  design  is  the  one  which  minimizes 

E  E  Var  (f*>  -  f<">)  - 

ges her 

among  all  designs,  under  a  0-way  ,  1-way  or  2-way  elimination  of  heterogeneity  model. 
Majumdar  (1986)  has  studied  this  problem,  and  has  identified  A-optimal  designs  in  various 
settings.  Here  we  give  an  example  of  an  A-optimal  design  for  one  way  elimination  of 
heterogeneity  and  an  example  of  a  design  which  is  optimal  for  each  of  0-,  1-  and  2-way 
elimination  of  heterogeneity. 


Example  1.  Suppose  4  test  treatments  are  to  be  compared  with  3  controls  in  30  blocks  of 
size  3  each,  under  1-way  elimination  of  heterogeneity  model.  Denoting  the  test  treatments 
by  A,B,C,D  and  the  controls  by  1,2,  and  3  the  A-optimal  design  is: 

112112112112111111 
233233233233  A  A  A  B  B  C 
AAABBBCCCDDDBCDCDD 


222222333333 

AAABBCAAABBC 

BCDCDDBCDCDD 

Example  2.  Suppose  8  test  treatments  are  to  be  compared  with  2  controls  in  a  12  x  12 

array  under  a  2-way  elimination  of  heterogeneity  model.  Denoting  the  test  treatments  by 

A,  B,C,D,E,F,G,H,  and  the  controls  by  1,2  ,  the  A-optimal  design  is  given  by: 
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This  design  is  model  robust  in  the  sense  of  being  A-optimal  for  0-  and  1-way  elimination 
of  heterogeneity  models  too. 

6.  Discussion. 

Cox  (1958,  p.  238)  advocated  augmenting  a  BIB  design  in  test  treatments  with  one 
or  more  replications  of  controls  in  each  block  as  a  means  of  getting  good  designs.  He 
neither  formally  mathematized  the  problem  nor  gave  any  justification  for  his  suggestion. 
However,  based  on  what  has  been  developed  during  the  past  several  years,  we  know  that 
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this  is  an  excellent  method  of  getting  efficient  designs  in  many  cases.  Pearce  (1960)  gave  a 
solution  for  A-optimal  designs  for  0-way  elimination  of  heterogeneity  model  and  proposed 
a  class  of  designs  for  comparing  test  treatments  with  a  control  and  gave  their  analysis 
for  the  1-way  elimination  of  heterogeneity  model.  Freeman  (1975)  studied  some  designs 
for  comparing  two  sets  of  treatments  for  the  2-way  elimination  of  heterogeneity  model. 
Pe§ek  (1974)  compared  a  BIB  design  with  an  augmented  BIB  design,  as  suggested  by  Cox 
(1958),  in  estimating  control  —  test  treatment  contrasts  and  noticed  that  the  latter  was 
more  efficient.  Das  (1958)  had  also  looked  at  augmented  BIB  designs. 

Bechhofer  and  Tamhane  (1981)  were  the  first  to  study  the  problem  of  obtaining  op¬ 
timal  block  designs.  However  their  optimality  consideration  was  neither  A-  nor  MV- 
optimality,  but  for  the  problem  of  obtaining  optimal  simultaneous  confidence  intervals 
under  a  1-way  elimination  of  heterogeneity  model.  Their  discoveries  led  to  the  concept  of 
BTIB  designs;  Notz  and  Tamhane  (1983)  studied  their  construction. 

Constantine  (1983)  showed  that  a  BIB  design  in  test  treatments  augmented  by  a 
replication  of  the  control  in  each  block  is  A-optimal  in  the  class  of  designs  with  exactly 
one  replication  of  the  control  in  each  block.  Jacroux  (1984)  showed  that  Constantine’s 
conclusion  remains  valid  even  when  the  BIB  designs  are  replaced  by  some  group  divisible 
designs. 

Majumdar  and  Notz  (1983)  gave  a  method  of  obtaining  A-optimal  designs  among  all 
designs  for  the  1-way  elimination  of  heterogeneity  model.  Hedayat  and  Majumdar  (1984) 
gave  an  algorithm  and  a  catalog  of  A-optimal  designs  and  studied  approximations.  Ture 
(1982)  also  studied  A-optimal  designs  and  their  approximations  and  construction.  Hedayat 
and  Majumdar  (1985a)  gave  families  of  A-optimal  designs.  Notz  (1985)  studied  optimal 
designs  for  the  2-way  elimination  of  heterogeneity  model.  Majumdar  (1986)  considered  the 
problem  of  finding  optimal  designs  for  comparing  the  test  treatments  with  two  or  more 
controls. 

Jacroux  (1984b,  1985)  gave  new  methods  for  obtaining  MV-optimal  designs  under 
1-way  elimination  of  heterogeneity  models,  gave  catalogs  and  studied  approximations. 
Jacroux  (1984c)  studied  optimal  designs  for  2-way  elimination  of  heterogeneity  models, 
utilizing  techniques  of  Hall  (1935)  and  Agarwal  (1966).  Hedayat  and  Majumdar  (1985b) 
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studied  designs  simultaneously  optimal  under  both  1-  and  2-way  elimination  of  hetero¬ 
geneity  models.  Jacroux  (1986)  generalized  the  Hedayat  and  Majumdar  (1984)  algorithm 
for  finding  A-optimal  designs.  Cheng,  Majumdar,  Stufken  and  Ture  (1986)  gave  new  fam¬ 
ilies  of  A-optimal  designs  and  some  approximations  for  1-way  elimination  of  heterogeneity 
models.  Stufken  (1986)  studied  A-optimal  designs  for  1-way  elimination  of  heterogeneity 
models,  gave  families  and  studied  approximations. 

There  are  many  other  design  settings  in  which  it  would  be  useful  to  identify  optimal 
designs  for  comparing  test  treatments  with  controls.  One  such  setting  is  that  of  repeated 
measurements  designs.  Some  aspects  of  optimality  and  construction  of  designs  in  this  area 
has  been  investigated  by  Pigeon  (1984). 

Giovagnoli  and  Wynn  (1985)  studied  A-optimality  of  designs  for  1-way  elimination  of 
heterogeneity  models  set  in  the  context  of  approximate  theory,  i.e.  with  an  infinite  number 
of  observations.  Spurrier  and  Edwards  (1986)  did  a  similar  study  for  optimal  designs  for 
finding  simultaneous  confidence  intervals. 

It  seems  appropriate  to  make  a  comment  on  randomization.  In  running  optimal 
designs  we  often  have  to  follow  a  well  structured  pattern.  This  does  not,  however,  mean 
that  there  will  be  no  room  for  randomization.  The  labelling  of  the  treatments,  experimental 
units  under  a  O-way  elimination  of  heterogeneity  model,  blocks  under  a  1-way  elimination 
of  heterogeneity  model  and  rows  and  columns  under  a  2-way  elimination  of  heterogeneity 
model  can  be  randomized. 
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